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We assay how inflationary models whose properties are dominated by the dynamics of a single scalar 
field are constrained by cosmic microwave background (CMB) data from the Wilkinson Microwave 
Anisotropy Probe (WMAP). We classify inflationary models in a plane defined by the horizon- flow 
parameters. Our approach differs from that of the WMAP collaboration in that we analyze only 
WMAP data and take the spectral shapes from slow-roll inflation rather than power-law param- 
eterizations of the spectra. The only other information we use is the measurement of h from the 
Ifubble Space Telescope (HST) Key Project. We find that the spectral index of primordial density 
perturbations lies in the la range 0.94 < Ws < 1.04 with no evidence of running. The ratio of the 
amplitudes of tensor and scalar perturbations is smaller than 0.61 and the inflationary scale is below 
2.8x10^® GeV, both at the 2cr C. L. No class of inflation or ekpyrotic/cyclic model is excluded. The 
A</!>^ potential is excluded at 3cr only if the number of e-folds is assumed to be less than 45. 



That there are multiple peaks in the CMB has re- 
cently been reinforced by data from the WMAP satel- 
lite [1-3]. This establishes that the curvature fluctua- 
tions which seed structure formation were generated at 
superhorizon scales. The inflationary paradigm [4] , which 
hinges on this very fact [5], is therefore vindicated. Other 
generic predictions of inflation, including approximate 
scale-invariance of the power spectra, flatness of the Uni- 
verse and adiabaticity and gaussianity of the density per- 
turbations are also fully consistent with WMAP data [6] . 

Given this supporting evidence for inflation, we adopt 
the full set of predictions of slow-roll inflation and obtain 
constraints on inflationary models imposed by WMAP's 
data. Studies of this type have been carried out in the 
past with less precise data [7] and with the simple pa- 
rameterization of the power spectra as power laws [8] . It 
has been emphasized that to extract precise information 
from data of WMAP's quality, high- accuracy predictions 
of the power spectra resulting from slow-roll inflation 
should be used [9]. There is a wealth of cosmological 
information in the CMB [10], and our approach is to em- 
ploy precise theoretical expectations in its extraction. 

Since the WMAP collaboration has considered what 
implications their data have for inflation [11], we describe 
at the outset the differing elements between our analy- 
sis and theirs. We elaborate on these differences later. 
WMAP include CBI [12], ACBAR [13], 2dFGRS [14] and 
Lyman-a power spectrum [15] data in addition to their 



own data. We restrict ourselves to WMAP data with 
a top- hat prior on the Hubble constant h {Hq — lOO/i 
km/s/Mpc), from the HST [16]. While we use the actual 
theoretical predictions for the primordial power spectra 
from single-field slow-roll inflation, WMAP parameter- 
ize the spectra with power-laws and a running spectral 
index. As a result, the spectral shapes we use are dif- 
ferent from those of WMAP and we directly fit to slow- 
roll parameters, while WMAP fit to derivative quantities. 
There are different virtues of the two approaches, and a 
comparison of the results obtained from them serve as a 
check of the robustness of the conclusions reached. 

Primordial spectra: 

The primordial scalar and tensor power spectra to 
0{lnk) are* [9,17] 

k k 
= A, (ao + ai In — + 02 In^ — ) , (1) 

Ph = At {bo + fei In ^ + 62 In^ A) , (2) 

where the pivot fc^ typifies scales probed by the CMB. 
The constants and bi are functions [9,18] of the 



*X is the intrinsic curvature perturbation and hij is the 
transverse traceless part of the metric tensor. 
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horizon-flow parameters, Cj of Ref. [19], that are defined 
by 
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Here, N is the number of e- folds since some moment ti 
during inflation, when the Hubble parameter was Hj. 

Note that to O(e^), the bi depend only on ei and £2, 
while the depend on ei, £2 and £3. We initially set 
£3 = 0, but wc then later demonstrate that the flt to 
the WMAP data is essentially not improved by includ- 
ing nonzero £3. The reason for not simply using the £3- 
indepcndcnt 0{e) expressions is that the 0{e'^) expres- 
sions are more accurate far from the pivot, and for a 
wider range of £1 and £2 [7,9]. It is uncertain that even 
high-precision data from the Planck satellite [20] can con- 
strain £3 to be small. Including £3 in our analysis would 
simply enlarge the allowed parameter space to include 
models which are not inflationary in the sense that the 
horizon-flow parameters are not small. 

The primary advantage of the horizon-flow parameters 
is that accurate predictions of the shapes and normaliza- 
tions of the power spectra can be made independent of 
parameters describing cosmic evolution. The horizon- 
flow parameters £1 and £2 are related to the usual slow- 
roll parameters [21] 
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via the flrst order relations [9] 

£l^£, £2^2(£-r?). (7) 

The normalizations of the spectra are given by 



A,= 



'PI 



16Hf 
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and the ratio of the amplitudes of the spectra at the pivot 

k — ^ IS 

Atbo 



R = = 16£l [1 + C£2 

+ (C^-Y+5)£,£2+(^-^ + l)£^] 



(9) 



where C = 7£;-|-ln2 — 2« —0.7296. Note that Uq, and 60 
are 0{1) and jao — 6o| is ©(£2). The spectral indices and 
their running can be expressed in terms of £1, £2 and £3: 

Tls = 1 - 2£i - £2 - 2el - (2C + 3)£i£2 - C£2£3 , (10) 
nt= -2£i-2£2-2(C+l)£i£2, (11) 

'^"^ =_2ei£2-£2£3, (12) 
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at = 



dink 

dnt 
dink 



= -2£l£2 . 



(13) 



We set £i = 0, i > 3. R, Us are rit are interdependent and 
the following consistency condition on single-fleld slow- 
roll inflation applies: 



R ~ -8nt 



(14) 



By our choice of formalism, we implicitly assume this 
condition to be satisfied. 

Note that the six inflationary parameters, As, At, ris, 
nt, as and at are determined by just three parameters. 
As, £1 and £2 in our analysis. In contrast, the WMAP 
collaboration parameterize the power-spectra with [11] 



(15) 



They eliminate rit as a free parameter by using the con- 
sistency condition Eq. (14). Thus, they have four free 
parameters, Ag, At, Us and as and set at = 0. 

A convenient classification of models based on the sep- 
arate regions in the R—Ug plane they populate, or equiva- 
lently relationships between the slow-roll parameters, was 
introduced in Ref. [22]. This classification becomes par- 
ticularly simple in the £2 — £1 plane, as discussed in the 
following section and shown in Fig. 3. 

Inflation models: 

In what follows wc use the common jargon, "red-tilt" 
for rig < 1 and "blue-tilt" for > 1. 
a) Canonical potentials of large-field models are the 
monomial potential, 



V{cP) = Vo{<l>/tir, p>2, 



(16) 



and the exponential potential V{(p) = Voe'^^^ of power- 
law inflation. They arc typical of chaotic inflation [23] 
and have V" > 0. The value of the scalar field falls 
0{Mpi) while the relevant perturbations are generated 
and thereby offer a glimpse into Planckian physics. In 
terms of the horizon-fiow parameters, large-field models 
satisfy 

< £2 < 4£i . (17) 

They have large R and predict a red tilt. 

b) Generic small-field potentials are of the form 

y(0) = Vo[l - (0//i)^] , P>2, (18) 

and are therefore characterized by V" < 0. The scalar 
field rolls from an unstable equilibrium at the origin to- 
wards a non-zero vacuum expectation value. Models re- 
lying on spontaneous symmetry breaking yield such po- 
tentials [24]. For small-field models 



£2 > 4£l 



(19) 



The tensor fraction is small and the scalar spectrum is 
red-tilted. 
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c) Potentials for hybrid inflation [25] are of the form 

V{cl>) = Vo[l + {Hl^r], P>2. (20) 

Hybrid inflation models involve multiple scalar fields. 
One of the fields, (f>, is the slowly rolling inflaton which 
does not carry most of the energy density {(f) < ^). An- 
other field which has a fixed value during the slow-roll 
of (f) provides Vq. When (f) falls below a critical value, 
the other field is destabilized and promptly ends infia- 
tion. As a result the value of (j) at the end of inflation 
is very model-dependent, which makes the number of e- 
folds correspondingly uncertain. Hybrid models can be 
treated as single-field models because the only role of the 
second field is to end inflation, and the slow-roll dynam- 
ics is dominated by a single field. These potentials arise 
in supersymmetric and supergravity models of inflation. 
Hybrid models have 

e2<0. (21) 

There is no robust prediction for R as can be seen from 
the above ei-independent inequality. However, a unique 
prediction of these models is that the spectrum is blue- 
tilted if 1 62 1 > 2ei. 

The line e2 =0 (i? = 8(1 — n,,)) implies e = r] which 
occurs for the exponential potential. Thus, power-law 
inflation marks the boundary between large-fleld and hy- 
brid models, 
c) Linear potentials, 

¥{<!>) = Voi<t>/l^), ¥{<!>) = Vo[l-{<t>/f^)], (22) 

define the boundary between large and small field models 
and lie at 

£2 = 4ei or 3R = 8(1 - n^) . (23) 

Since ei is a constant for such potentials, inflation ends 
only with the help of an auxiliary field or some other 
physics. 

To avoid overstating the comprehensiveness of this 

classification of potentials, we list a few potentials of dif- 
ferent forms [26], which however do fit into the large- 
field, small-field, hybrid or linear categories according to 
the relationships between the horizon-flow parameters. 
V{<j>) = yo[l±ln(</'//x)], Voil-e-^/'^] and Vo[l±{^/ h)-p] 
are hybrid in the sense that an auxiliary field is needed 
to end inflation, but lead to £2 > 4ei and therefore lie in 
the small-field region of the parameter space. Similarly, 
power-law inflation does not end without a hybrid mech- 
anism, but lies in the large-field region. Finally, let us 
note that the simplest models of the ekpyrotic/cyclic [27] 
variety are small- field [28]; however, the prediction for rig 
in these models is controversial [29]. 

Analysis: 



We compute the TT and TE power spectra ST^ = 
l{l + 1)Ci/2tt, using the Code for Anisotropies in the 
Microwave Background or CAME [30] (which is a paral- 
lelized version of CMBFAST [31]) and a supporting mod- 
ule that calculates the inflationary predictions for the pri- 
mordial scalar and tensor power spectra [32] . We assume 
the Universe to be flat, in accord with the predictions of 
inflation, and that the neutrino contribution to the mat- 
ter budget is negligible. The dark energy is assumed to 
be a cosmological constant A. We calculate the angular 
power spectrum on a grid consisting of the parameters 
specifying the primordial spectra^, ei, e2 and ^4^, and 
those specifying cosmic cvolution-t , the Hubble constant 
h, the reionization optical depth r, the baryon density 
Ub = flsh^ and the total matter density ojm = Om/i^ 
(which is comprised of baryons and cold dark matter). 
We choose h rather than f^A because it is directly con- 
strained by the HST [16]. We do not include priors 
from supernova [33], gravitational lensing [34] or large 
scale structure [35] data or nucleosynthesis constraints 
on Wb [36] although these would somewhat sharpen the 
cosmological parameter determinations. 

We employ the following top-hat grid: 

• 0.0001 < ei < 0.048 in steps of size 0.004 

• -0.18 < £2 < 0.14 in steps of size 0.02 

• 0.64 <h< 0.80 in steps of size 0.02 

• r =0, 0.05, 0.1, 0.125, 0.15, 0.16, 0.17, 0.18, 0.19, 

0.2, 0.25, 0.3 

• 0.018 <ujt< 0.028 in steps of size 0.001 

• 0.06 <(Jm < 0.22 in steps of size 0.02 

• As is a continuous parameter. 

The range of values chosen for h correspond to the HST 
measurement, h = 0.72 ± 0.08 [16]. This serves to break 
the degeneracy between Om and Oa without the need for 
supernova data [37]. We place the pivot at k^, = 0.007 
Mpc~^. The primordial spectra are evaluated to 0{\n^ k) 
(see Eqs. 1-2) with £3 = 0. 

The WMAP data are in the form of 899 measurements 
of the TT power spectrum from I = 2 io I = 900 [2] and 
449 data points of the TE power spectrum [3]. We com- 
pute the likelihood of each model of our grid using Ver- 
sion 1.1 of the code provided by the collaboration [38]. 
The WMAP code computes the full covariance matrix 
under the assumption that the off-diagonal terms are sub- 
dominant. This approximation breaks down for unreal- 
istically small amplitudes. We include the off-diagonal 



^We call these inflationary parameters. 
*We call these cosmological parameters. 
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terms only when the first peak occurs between I = 100 
and I = 400 with a height above 5000 fiK^ [38]. The 
restriction on the peak location may at first appear ir- 
relevant, but for very large tensor amplitudes, the max- 
imum height of the scalar spectrum shifts to very small 
I. When the height of the first peak is below 5000 ^iK^ 
(which is many standard deviations away from the data) , 
we only use the diagonal terms of the covariance matrix 
to compute the likelihood. To obtain single-parameter 
constraints, we plot the relative likelihood e'^'^™*""-^ 
for each parameter after maximizing over all the others. 
The x-a range is obtained for likelihoods above 
For 2-dimensional constraints, the la, 2a and 3cr regions 
are defined by Ax^ = 2.3, 6.17 and 11.83, respectively, 
for two degrees of freedom. 

Results: 

Although our primary focus is to obtain constraints 
on the inflationary models, we display constraints on the 
cosmological parameters as a consistency check with the 
WMAP collaboration. This is pertinent because we have 
used a form for the primordial spectra (Eqs. 1-2) that is 
specific to single-field slow-roll inflation rather than the 
standard power-law form. 

The best-fit parameters are ei = 0.016, £2 = —0.02, 
u)b = 0.023, ujm = 0.12, T = 0.125 and h = 0.78 with 
normalization ^ 22.0 x 10"!° with = 1428.88 for 
1341 degrees of freedom^. The results from Fig. 1 are 
in good agreement with those obtained by the WMAP 
collaboration [6]. This indicates that the choice of the 
spectral shapes does not matter at the present precision 
of the WMAP data. 

Likelihood plots for the inflationary parameters are 
shown in Fig. 2. We do not show the result for rif because 
we have imposed the consistency condition, Eq. (14). The 
Icr confidence limits on various parameters are provided 
in Table. I. We see that the spectra are consistent with 
scale-invariance and with a small tensor contribution; the 
best-fit scale-invariant spectrum with no tensor contribu- 
tion has = 1430.61. Also, running of the spectral in- 
dices is insignificant. (In our framework OLg^t are required 
to be consistent with inflationary predictions as dictated 
by Eqs. (12) and (13); they are not free parameters). 

WMAP has provided important information about R 
and the energy scale of inflation. 



R < 0.61 (2c7 limit) , 



(24) 



and 



^To check the validity of setting = 0, t > 3, we enlarged 
the grid to include £3 = ±0.1 and found the minimum 
values to be ^{tz = -0.1) = 1428.80 and x^(e3 = +0.1) = 
1429.05. The small changes in the x^ values confirm that it 
is not necessary to include es in the analysis. 
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FIG. 1. Relative likelihood plots for several cosmological 
parameters and Ag. to is the age of the Universe in Gyr. We 
do not show the plot for h because it is not constrained by 
the fit; see Table. I. 
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or equivalently, 



= y^neiAs < 1.48 X 10"^ {2a limit) , (25) 
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Mpi < 2.8 X 10^'' GeV {2a limit) . 

(26) 



Since V/ > 10^^ GeV is consistent with data, it is still 
possible to detect inflationary gravity waves by measur- 
ing the curl component of CMB polarization [39]. 

Implications for models: 

The allowed regions of ei and 62 (and equivalently Hg 
and R) are shown in Fig. 3. The different classes of in- 
flationary models populate distinct regions of the £2 — ei 
plane, as discussed above. The consistency of the models 
with the data can be judged from the figure. Even with 
the high quality of the WMAP data, no class of models 
is excluded. As long as ei <C 1, £2 = is consistent with 
data, this will remain the case. Moreover, the allowed 
range for (see Table I) is consistent with the predic- 
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FIG. 2. Relative likelihood plots for some inflationary pa- 
rameters. 



tions of a wide spectrum of inflationary models, and so 
does not help in discriminating between them. 

We now place some constraints on large-field and 
small-field models whose predictions do not involve too 
much freedom. 

For the monomial potentials (p > 2) of large-field mod- 
els, 



V 



N ~ 



V 



(27) 

(28) 



where N is the number of c-folds of inflation from the 
time that scales probed by the CMB leave the horizon 
until the end of inflation**. At least about 40 e- folds are 
needed for the Universe to be flat and homogeneous, and 
typically the largest value is 70 [40]. Since ei <C 1, £2 = 
is allowed, p cannot be constrained independently of N . 

The 3cr exclusion of the potential in Ref. [11] was 
based on an analysis of WMAP data in combination with 
higher I CMB data and large scale structure data, as- 
suming = 50 (for which = 0.94 and R = 0,32), 
Their results are shown in the second row of their Fig. 4. 
Note that the point ris = 0.94, R = 0.32 lies inside the 
3(T region of our Fig. 3 and is therefore not excluded by 
WMAP data alone. If instead N is 60, then = 0.95 



**We are reverting to the conventional definition in which 
the number of c-folds is counted backward in time; in the 
definition of the horizon-fiow parameters, e-folds are counted 
forward in time. 
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TABLE I. The la limits on the inflationary and cosmo- 
logical parameters. Quantities below the line axe not directly 
constrained by the data but derived from those above the line. 



and R = 0.27; this point lies in the 95% C. L. allowed 
region of the second row of Fig. 4 of Ref. [11] and within 
the 2(7 region of Fig. 3. 

If ei ^ (e2 > 0), a lower bound (upper bound) p/4 > 

^niin^^max (^^4 ^ ^max ^ ^rntn-^ cusucs; p is determined if 

both conditions on the e's are simultaneously satisfied. 
Since we expect p <C 4A'', 



N ^ 



^max 
^2 



12.5 {2a limit) , 



and 



p < ANeT"'' = 0.15N {2a limit) . 
For small-field models with p > 3 [26] , 

2 p~l 



Np-2 



(29) 



(30) 



(31) 



The least stringent bound on N occurs in the limit p 
oo. 



N > 



^2 



25 {2a limit) . 



For the small- field quadratic potential {p = 2), 



ei -C €2 = 



27r 



We find the 2a bound 



7nax 
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1AM PI . 



(32) 



(33) 



(34) 



Similar constraints can be placed on other models, but 
unfortunately, they are not very enlightening. 

Conclusions: 
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FIG. 3. la, 2a and 3a allowed regions in the £2 — ei and 
R — Us planes. We have plotted the predictions for the Xif)^ 
potential with the number of e-folds A'' = 40, 50, 60 and 70. 
The prediction approaches ei = £2 = (us = 1, = 0) as 

^ oo. 



WMAP has provided compelling evidence for the in- 
flationary paradigm. We have adopted the explicit pre- 
dictions of single-field slow-roll inflation for the shapes of 
the power-spectra to analyze WMAP data. The fact that 
our parameter determinations are consistent with those 
obtained with the standard power-law parameterization 
by the WMAP collaboration provides further evidence 
for slow-roll inflation. Since exact scale-invariance and 
a negligible tensor contribution to the density perturba- 
tions are adequate to describe the data, it is not presently 
possible to exclude classes of inflationary models. 

We have shown how different classes of inflationary 
models can be distinguished in the €2 — ei plane of 
the horizon-flow parameters. If and when the horizon- 
flow parameters Ci or/and £2 are determined to be non- 
zero, large numbers of inflationary models will be ruled 
out. For that, we await even higher precision data from 



WMAP and eventually from Planck. 
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